BCS-BEC Crossover in Symmetric Nuclear Matter at Finite Temperature: Pairing Fluctuation and 

Pseudogap 
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By adopting a T-matrix-based method within the GqG approximation for the pair susceptibility, we studied 
the effects of pairing fluctuation on the BCS-BEC crossover in symmetric nuclear matter. The pairing fluctuation 
induces a pseudogap in the excitation spectrum of nucleon in both superfluid and normal phases. The critical 
temperature of the superfluid transition was calculated. It differs from the BCS result remarkably when density 
is low. We also computed the specific heat, which shows a nearly ideal BEC-type temperature dependence at 
low density, but a BCS-type behavior at high density. This qualitative change of the temperature dependence of 
specific heat may serve as a thermodynamic signal for the BCS-BEC crossover. 

PACS numbers: 21.65.-f, 21.30.Fe, 26.60.-c, 74.20.-z 



I. INTRODUCTION 

One of the most common properties of attractive fermion 
many-body system is the arising of superfluid state at low 
temperature. Depending on the strength of the attractive in- 
teraction between two fermions, however, the physical con- 
tents of the superfluid state could be distinguishably altered. 
When the interaction is weak, the system can be well de- 
scribed by the Bardeen-Cooper-Schrieffer (BCS) theory. In 
this case, the superfluidity is due to the condensate of loosely 
correlated Cooper pairs and superfluid gap is much smaller 
than the Fermi energy. When the interaction becomes suffi- 
ciently strong, the two-fermion bound state could form, which 
may behave like a boson. In this situation, the superfluidity 
is due to the Bose-Einstein condensation (BEC) of the tightly 
bound two-fermion state and the superfluid gap could be much 
larger than the Fermi energy. Although the BCS and BEC 
limits have quite different physics, it was found that there is 
no true phase transition (traditional symmetry breaking) hap- 
pening in between. The transition from BCS state to BEC 
state is smooth and is often called the BCS-BEC crossover [1- 
4]. Such BCS-BEC crossover was recently realized in cold 
atomic experiments (see Refs. [5, 6] and references therein). 

It was well known that cold nuclear matter can be in a su- 
perfluid state, which plays a crucial role in a variety of nuclear 
many-body problems, from neutron stars, low-energy heavy- 
ion collisions, to finite nuclei. It was argued that the BCS- 
BEC crossover should also occur in nuclear matter where the 
BCS state of neutron-proton (tip) Cooper pairs at high density 
undergoes a smooth transition into BEC state of deuterons at 
low density [7-14]. At the same time, the chemical potential 
changes its sign at a certain density and finally approaches 
one-half of the deuteron binding energy at the low-density 
limit. Recently, a similar situation was also studied for the 
neutron-neutron (mi) pairs in the 1 S'o channel [14-25]. It was 
found that in certain (low) density regions the nn pairs can 
be strongly correlated. However, no assured BEC state was 
found for nn pairs. 

So far, most of the investigations of nuclear BCS-BEC 
crossover in the literatures focused on the ground state 



crossover described by BCS theory. Although the BCS the- 
ory succeeds in describing the BCS-BEC crossover at zero 
temperature, it, as a mean-field theory, is not sufficient to de- 
scribe low-density nuclear matter at finite temperature where 
the pairing fluctuation is substantial due to the strongly cor- 
relating nature of the system. Actually, as a consequence of 
the strong correlation, the low-density nuclear matter exhibits 
"pseudogap" phenomena above the critical temperature T c of 
superfluid transition and has an exotic normal state that is 
different from the Fermi liquid normal state associated with 
BCS theory [26-28]. Similar situations were also found in 
other strongly correlated systems, such as high T c supercon- 
ductors [29-34] and cold atomic Fermi gases under Feshbach 
resonances [33, 35-40]. To include the pairing fluctuation ef- 
fects and investigate the pseudogap phenomena, we will adopt 
a T-matrix formalism based on a Go G approximation for the 
pair susceptibility which was first introduced by the Chicago 
group [34^-0]. This formalism generalizes the early works 
of Kadanoff and Martin [41] and Patton [42], and can be con- 
sidered as a natural extension of the BCS theory since they 
share the same ground state. Moreover, this formalism allows 
quasi-analytic calculations and gives a simple physical inter- 
pretation of the pseudogap phase. It clearly shows that the 
pseudogap is due to the incoherent pairing fluctuation. 

Our focus will be put on the np pair in symmetric nuclear 
matter (mainly in the low-density region), since the interac- 
tion in this case is more attractive than in the nn or pp chan- 
nels and it provides a very good playground for the BCS-BEC 
crossover. We will extend the early studies [7-14] to include 
the pairing fluctuation effects and determine the magnitude of 
the pseudogap. Furthermore, the transition temperature for 
the onset of the superfluid and the thermodynamic properties 
will be also be a concern. Such a study will be helpful to un- 
derstand the strongly coupling nature of low-density nuclear 
matter and may give useful information on the physics of the 
surface of nuclei, expanding nuclear matter from heavy-ion 
collisions, collapsing stars, etc. 

The article is organized as follows. We give a brief sum- 
mary of the effective nucleon-nucleon potential in Sec. II. In 
Sec. Ill, we give a detailed theoretical scheme of how the T- 
matrix-based formalism works at finite temperature. The nu- 
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merical results are presented in Sec. IV. We summarize our 
results in Sec. V. Throughout this article, we use natural units 

h = ks = c = 1, 



alized mainly in the spin-triplet s-wave channel, so let us con- 
sider the following Lagrangian describing neutron and pro- 
ton interaction via two-body attractive forces in 3 Si, S z = 
channel, 



II. EFFECTIVE NUCLEON-NUCLEON POTENTIAL 

The aim of this article is not to determine the precise values 
of the pairing gap, the critical temperature, etc., but rather to 
perform a qualitative (or semi-quantitative) study of the ef- 
fects of pairing fluctuation on BCS-BEC crossover. In or- 
der to highlight the essential physics, we will adopt a sim- 
ple density dependent contact interaction (DDCI) developed 
in Refs. [43, 44]. The potential is of the form 



V{x -x')=vo\l-n 



Pq 



5(x-x'), (2.1) 



where t>n,f7, 7 are three adjustable parameters, p(x) = 
p n (x) + p p (x) is the nuclear density and po = 0.17 fin 
is the normal nuclear density. Taking suitable values of the 
parameters, one can reproduce the pairing gap A(fcp) as a 
function of Fermi momentum kp = (3tt 2 p/2) 1 ^ 3 in the chan- 
nels L = 0,1 = l,I z = ±1,5 = and L = 0,1 = 
0, S = l,S z = calculated from realistic nucleon-nucleon 
potentials [43, 44], where L is orbital angular momentum, 
I denotes isospin, and S is spin. According to Garrido 
et al. [43, 44], we will choose in the following numerical cal- 
culation r) = 0, v = -530McVfm 3 in the 1 = 0, 3 5i (np 
pairing) channel and a energy cutoff e c = 60 McV to regu- 
larize the integration. With these parameters one must use a 
density-dependent effective nucleon mass m(p) correspond- 
ing to the Gogny interaction[43, 44], 



m(p) 



m 



L i + ^hlL j^[ Wc + 2(B C H c ) - AM C 



cosh: 



sinh : 



(2.2) 



where x c = k 2 F p 2 /2, mo = 939 MeV is the bare mass of 
nucleon, and p c , W c , B c , H c , M c are parameters correspond- 
ing to the Gogny force Dl[45, 46], their values are listed in 
Table I. 



TABLE I: Parameters in the effective mass of nucleon (2.2) corre- 
sponding to the Gogny interaction Dl[45, 46]. 



c 


jii c [fm] 


!V c [MeV] 


B c [MeV] 


H c [MeV] 


M c [MeV] 


1 


0.7 


-402.4 


-100.0 


-496.2 


-23.56 


2 


1.2 


-21.30 


-11.77 


37.27 


-68.81 



m. T-MATRIX-BASED FORMALISM 

We consider the nuclear matter as an infinite system of in- 
teracting fermions. In the low-density region, np pairing is re- 



c = $*,° ( - d - + 2^ + M ) ^ 



i—,np a— "f ,4, 

(3.1) 

where g = — Vq > is the coupling strength in np chan- 
nel and t = it is the imaginary time. Introducing aux- 
iliary fields A = (g/2) (</> n tA4 ~ ^pf^ni) and A = 
(g/2) (tppi^n-r ~ i>nl4>rt), we can recast Eq. (3.1) as 

C = V V 1pt,a \-d T + — + A« ll>i,a ~ -AA 

^-^ \ 2m ] q 

i=,np cr=t,l 

+ A {ippllpnf ~ lpnli>pt) + {4>nt4>pl ~ VVfVw) A, 



-AA, 
9 



(3.2) 



where we have introduced the Nambu-Gorkov spinor ^ = 

("0nt , ^pI , ippt i ^ni) T , and 



:-i - (Si 1 



s- L = 



S^ 1 



(3.3) 



with 



c _i _ f-d T + V 2 /{2m) + p A 

"1=1 & a Y72 



-d r + V 2 /(2m) + /i 



-d T - V7(2m) - fi 



-d T - V 2 /(2m) - n 



(3.4) 



(3.5) 



It is seen that 52 differs from S\ only by minus signs in front 
of A and A. To make our formulae more compact, in the 
following discussions we will treat S\ only and neglect the 
subscript 1 without confusion. 

In the rest of this section, following the works of the 
Chicago group, we will introduce the basic method of the T 
matrix. This T matrix is defined as an infinite series of ladder- 
diagrams in particle-particle channel (rather than particle-hole 
channel) by constructing the ladder by one free nucleon prop- 
agator and one full nucleon propagator. Then, as usual, the T 
matrix enters the nucleon self-energy in place of the bare in- 
teraction vertex. The coupled T-matrix equation and the self- 
energy equation (as well as number density equation) should 
be solved self-consistently. One can view this approach as the 
simplest generalization of the BCS scheme, which formally 
can also be cast in a T-matrix formalism. Let us discuss this 
point in the following subsection. 



3 



A. BCS theory 

The BCS theory is based on mean-field approximation to 
the (anomalous) self-energy, i.e., A and A are chosen as their 
mean-field values A = A s f and A = A s f (without loss of 
generality, we put A s f and A s f to be constants and A s f = 
A s f), which are regarded as order parameters for superfluid 
phase transition. We start with the Nambu-Gorkov formalism 
in momentum space, 



S mf^)-{ Asf -Q-\-K) 



(3.6) 



where K = (iu! n , k) and ioj n = i(2n + 1)itT is the fermion 
Matsubara frequency. Q^^K) = iu> n — £k is the inverse of 
free nucleon propagator, £k = k 2 / (2m) — /x is the dispersion 
relation of free nucleon. From Eq. (3.6) one gets, 



f9mi(K) T mi (K) 



\T mi {K) -G mi {-K) 
where F m f(K) is the anomalous propagator, 



(3.7) 



T ni {(K) = A sf g mf (K)g (-K) 
-A sf 

El' 



{iu n ) 2 

and Gmf(K) is the mean-field single nucleon propagator, 



(3.8) 



(iw„) 2 - 



(3.9) 



w ith the m ean-field dispersion relation of nucleon 
+ A 2 f and the mean-field self-energy 



Smf(A") = -A 2 sl G Q {-K). 
The coupled gap and density equations read 



(3.10) 
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K 



gA sf 
V 



2n F (£ k )] , 



_ V p tr > u 



-E 

v 4^ 



1 - -g (1 - 2n F (£ k )) 



where n p (x) — 1/ [cxp (f3x) + 1] is the Fermi-Dirac function 
and e"' IUJn with r) —> is a convergence factor for Matsubara 
summation. The prefactor 2 on the right-hand side of density 
equation counts the degeneracy of Si and S2- 

In BCS theory, np pairs enter into the problem below T c , 
but only through their condensates at zero momentum. By 
rewriting the mean-field self-energy S m f (K) in a manner of 

£mf(*0 = ^ $>mf(Q)So(Q - K), 

P Q 



we are aware of that these condensed pairs can be associated 
with a T matrix in the following form 



t ml {Q) = -A 2 sl S(Q). 



(3.13) 



with Q = (qo, q), qo = iuj v = ilvitT, v 6 Z being the boson 
Matsubara frequency and 5(Q) = f35^ y oS^' (q). Furthermore, 
if we define the mean-field pair susceptibility as 

xMQ) = jjyY,GMK)Go(Q-K), (3.14) 



K 



we can re-write the gap equation in superfluid phase as 

l-5Xmf(0) = 0, T<T C . (3.15) 

This suggests that one can consider the uncondensed pair 
propagator or T matrix to be of the form 



1 - 9Xmt(Q) 



(3.16) 



and then the gap equation is given by t~ ir (0) = 0. 

It is well known that the critical temperature T c in the BCS 
theory is related to the appearance of a singularity in a T ma- 
trix in the form of Eq. (3.16) but with A S £ = 0. This is the 
so-called Thouless criterion for T c [47]. But the meaning of 
Eq. (3.15) is more general as stressed by Kadanoff and Mar- 
tin [41]. It states that under a asymmetric choice of x> the 
gap equation is equivalent to the requirement that the T ma- 
trix associated with uncondensed pair remains singular at zero 
momentum and energy for all temperatures below T c . 

Although the construction of the uncondensed pair propa- 
gator (3.16) in BCS scheme is quite natural, the uncondensed 
pair has no feedback to the nucleon self-energy (3.12). When 
the coupling is weak, such a feedback is not important, but if 
the system is strongly coupled, this feedback will be signifi- 
cant. The simplest way to include the feedback effects is to 
replace i m f in Eq. (3.12) by t m { + t pa ir- But to make such an 
inclusion self-consistent, i pa j r should be somewhat modified 
which we discuss in next subsection. 



B. G G formalism at T < T c 

(3.11) Physically, the BCS theory involves the contribution to nu- 
cleon self-energy below T c only from those condensed pairs, 
i.e., the q = Cooper pairs. This is justified only at weak- 
coupling region. Generally, in superfluid phase, the self- 
energy consists of two distinctive contributions, one from the 
superfluid condensate, and the other from thermal pair exci- 
tations. Correspondingly, it is natural to decompose the self- 
energy into two additive terms 

W = Z)*(G)So(Q - K) = S mf (JC) + Spg(-ftT), 

(3.12) Q 



(3.17) 
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with the T matrix accordingly given by 

t(Q) = t m{ (Q) + t pg (Q), 

tpg(Q) = 



(3.18) 



1 - 9X(Q) ' 

where the subscript pg indicates that this term will lead to 
the pseudogap in nucleon dispersion relation as will become 
clear soon. See Fig. 1 for the Feynman diagrams for t pg (Q) 
and E(A'). Comparing with the BCS scheme, t m {(Q) in 
Eq. (3.12) is replaced by t(Q), and now contains the 

feedback of uncondensed pairs. The pair susceptibility x(Q), 
as inspired by Eq. (3.14), is chosen to be the following asym- 
metric GqG form, 

= ^E^^WQ-*)- (3 - 19) 

In spirit of Kadanoff and Martin, we now propose the su- 
perfluid instability condition or gap equation as [extension of 
Eq. (3.15)] 



l-ffX(0)=0, T<T C . 



(3.20) 



We stress here that this condition has quite clear physical 
meaning in BEC regime. The dispersion relation of the bound 
pair is given by t~ 1 (Q) = 0, hence £ -1 (0) oc fib with fhb 
the effective chemical potential of the pairs. Then the BEC 
condition requires fib = for all T < T c . 

The gap equation (3.20) tells us that t pg (Q) is highly 
peaked around Q = 0, so we can approximate E pg as 

E PS (K) ~ -A 2 ps g (-K), T<T C , (3.21) 
where we have defined the pseudogap parameter via 

i-£t pg (Q). (3.22) 



PV 



The total self-energy now is 

E(JT) = -A 2 Q (-K), 



(3.23) 



with A 2 = A 2 f 



A 2 . 

pg 



It is clear that A D<r also contributes 



to the energy gap in quasi-nucleon excitation. Physically, 
the pseudogap A pg below T c can be interpreted as an extra 
contribution to the excitation gap of nucleon quasi-particle: 
an additional energy is needed to overcome the residual at- 
traction between nucleons in a thermal excited pair to pro- 
duce fermion-like quasi-particles. One should note that the 
Apg is associated with the fluctuation of the pairs A 2 



pg 



(A 2 ) - (A) 2 [34, 37], hence it does not lead to superfluid 
(symmetry breaking). 

With Eq. (3.21), the pair susceptibility reads, 



x(Q) 



2Eu 



2E k £ k + £ q - k -q - *0+ 

-Ek-Ck n F (E k ) - n F (£ q _ k ) 
£k - Cq-k + qo + i0 + 
sE k + £ k n F (-sE k ) - n F (£ q _ k ) 



1 E 



2sE k 



sE k + £ q _ k - qo 



(3.24) 







+ 


< 

— *- 





+ 



J l)g 



J nif 



FIG. 1 : Feynman diagrams for the T matrix of non-condensed pairs 
and the nucleon self-energy in the GoG formalism. 



with i?k = \/£ k + A 2 . The number equation remains un- 
changed except the replacement of A m f — » A. 

Furthermore, the gap equation (3.20) suggests that we can 
make the following pole approximation to the pair propagator 
or T matrix t ps (Q) as 



t PS (Q) 



z- 



Qo - q 2 /(2m b ) : 



(3.25) 



where the residue Z 1 and effective "boson" mass are given 
by 



(3.26) 



We stress here that in general, the expansion of t pg (Q) should 
also contain a term oc q 2 ,, but such term does not bring quali- 
tative change to the crossover physics [38], hence we neglect 
itinEq. (3.25). 

A straightforward calculation gives 



z = 


dx 








dqo 


Q=0 


z 


ld 2 X 






3<9q 2 


Q=0 



s n F (E k ) - n F (s£ k ) 



V 
1 

A 2 



k s=± 
P 
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2£ k 



E k - s£ k 



^E ni? ^ k ) 



(3.27) 



and 

Z 
mb 



^EE 

k s 



2sE k 



1 n F (-sE k ) - n F (£ k ) 



sEy 



2k 2 /n F (- S £ k )-n F (£ k ) 



3m 2 
Z 1 
m 



(sE k + Ck) 2 
2k 2 



V ^ 3m 2 A 2 

k 



^ sE k + £ k 
1 ^ k 2 



V ^ 3m 2 E k A i 

k 



x {(^k + £k)[l - 2n F (^ k )] - 2SkCk[l - 2njr(&)]} . 

(3.28) 

The expression in the square bracket of right-hand-side of 
Eq. (3.27) is nothing but one-half the density of the pairs Pb/2, 
we then have pb = 2ZA 2 . 
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Substituting Eq. (3.25) into Eq. (3.22) leads to 
A ™ = ^E n *fc a /(2m 6 )] 



I ( Tm b\ 

Z V 2tt J 



3/2 



(3.29) 



where tib(x) = l/[cxp (f3x) — 1] is the Bose-Einstein func- 
tion and a vacuum term was regularized out. It should be 
stressed that at zero temperature Ap g = 0, hence the GoG 
scheme yields the BCS ground state. One should note that 
A ps = Pb incondcnscd /2Z, and hence = ^condensed / 2 £. 

Now, Eq. (3.20), Eq. (3.29), as well as number equation are 
coupled to determine the total excitation gap A, the pseudo- 
gap Apg and the nucleon chemical potential fj, at given density 
and temperature below T c . In short, they are 



1 - fEin- 2 ^)] 



-E 



^(l-2n F (£ k )) 



3/2 



A 



I ( Tm A 

Z\2tt) ^ 



(3.30) 



C. GoG formalism above T c 

Above T c , Eq. (3.20) does not apply, hence Eq. (3.21) no 
longer holds. To proceed, we extend our more precise T <T C 
equations to T > T c in the simplest fashion. We will con- 
tinue to use Eq. (3.23) to parameterize the self-energy but with 
A = Apg, and ignore the finite lifetime effect associated with 
normal state pairs. It was shown that this is still a good ap- 
proximation when temperature is not very high [36, 38, 40]. 
The T matrix i pg (Q) at small Q can be approximated now as 



where Li n (z) is the polylogarithm function. Then Eq. (3.32), 
Eq. (3.33) and the number equation which remains unchanged 
determine A pg , fi and fib- 
In summary, at T > T c , the order parameter is zero, and 
A = Apg. The closed set of equations determining A, \i and 



M6 is 



Zfi b 



l r l- 2n F {E u ) 



2£ k 



P = 77 



V ^ 

k 



^ (1 - 2n F (E k )) 



A^ 



lf^V /2 Li3( 
Z \ 2tt / 3 V 



(3.34) 



D. Thermodynamics 

The thermodynamics of the matter are governed by the ther- 
modynamic potential, which reads 



(3.35) 



where Clf and Of, are the contributions from nucleons and 
thermal excited pairs, 



n 



2AM0)-f E 



V 



Ek — £k 



ln(l 



n h = 



pv 



E ln (! 



(3.36) 
(3.37) 



Other thermodynamic quantities can be derived from f2, e.g., 
the entropy density is given by s = —dfl /dT, and the specific 
heat cy is given by cy = Tds/dT. 



IV. NUMERICAL RESULTS 



WQ) 



q - Sl q 



(3.31) 



where J7 q = q 2 /(2mb) — /if,. Since there is no condensation 
in normal state, the effective pair chemical potential \x\, is no 
longer zero, instead, it should be calculated from 

Zfi h = t-\o) = --+x(0) 
9 



1 1^1- 2n F (E u ) 
9 



2E k 



This is used as the modified gap equation. Similarly, above T c 
the pseudogap A pg is determined by 



A p g = ^E**^ 



l_ ( Tm b \ 
Z \~2^ 



\ 3/2 

J Lis (e w ' /T ) , (3.33) 



We discuss now the results obtained by numerically solving 
Eq. (3.30) for T < T c and Eq. (3.34) for T > T c . We will 
mainly focus on the intermediate (the crossover region, see 
Fig. 2) and low density regions, since the high density region 
is proven to be well understood in BCS theory. We begin with 
the results concerning the critical temperature for superfluid 
transition in the BCS-BEC crossover. 



A. BCS-BEC crossover and critical temperature 

At zero temperature, the GoG formalism reproduces the 
usual BCS theory. In order to have a quantitative examina- 
tion of the BCS-BEC crossover, it is convenient to define the 
condensed rip Cooper pair wave function at zero temperature, 



V(r) = G^GSKVx^x 
d 3 k 



■t)\BCS) 



c 



(27T) 



M)e 



?k r 



(4.1) 
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where a' nc (at pa ) is the creation operator of neutron (proton) 
with spin a and ^>(k) are the anomalous density distribution 
function 



V(k) = (BCS\al t (k)al(-k)\BCS) 
A 



(4.2) 



After substituting Eq. (4.2) into the number and gap equations, 
we get the following Schrodinger-like equation, 



-V(k)- 5 (l-2n k ) 



^(k') = 2^(k). 



(4.3) 



In the limit of vanishing density, —> 0, this equation goes 
over into the Schrodinger equation for the np bound states 
(the deuterons) in the center-of-mass frame, and the chem- 
ical potential 2p then plays the role of the binding energy. 
Hence, one expects that at sufficiently low density and low 
temperature, the symmetric nuclear matter should be in the 
BEC phase. 

To have more quantitative description of the BCS-BEC 
crossover, we define other characteristic quantities: the mean- 
square-root size of the np pair, 



J d 3 xx 2 |V>(x)|' 

f Cpx|V(x)| 2 



(4.4) 



and the s-wave scattering length a that relates the coupling 
constant g to the low-energy limit of the two-body T matrix 
of np scattering in vacuum, 



m 
Ana 



1 

9 



d 3 k 



(2tt) 3 k 2 ' 



(4.5) 



In BCS region, £ is expected to be larger than the average dis- 
tance between neutron and proton d = (p/2) -1 / 3 and at the 
same time the scattering length a should be negative to ensure 
that the interaction between neutron and proton is attractive. 
In the BEC region, however, £/d should be very small reflect- 
ing the compactness of the pair, and the scattering length will 
be positive to guarantee the appearance of two-body bound 
state. 

In Fig. 2, we show the nucleon chemical potential over 
Fermi energy ratio p/ep, the scaled mean-square-root size 
£/d, and the inverse scattering length l/(kpa) as functions of 
nuclear density. Although £ itself is not a monotonous func- 
tion of p, the scaled one goes down monotonously with de- 
creasing density and finally approaches zero at zero density. 
The right vertical line around p/po ~ 0.5 indicates the posi- 
tion of £ = d which can be used to separate the BCS (weak 
coupling) region from the crossover (intermediately strong 
coupling) region. The chemical potential roughly equals the 
Fermi energy at BCS region, but it drops down with decreas- 
ing density and becomes negative below p/po ~ 0.002. The 
position where p changes sign can be regarded as the bound- 
ary between BEC (strong coupling) region in which p is neg- 
ative and other region with positive p. The third panel shows 
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FIG. 2: (Color online) The nucleon chemical potential over Fermi 
energy ratio p/sf, the scaled pair size and the scaled scattering 
length of np collision l/(kfO) as functions of density. The right 
and left vertical lines that separate the BCS, BEC, and crossover re- 
gions are, respectively, determined by the conditions £/d = 1, and 
p/sF = 0. The red vertical line in crossover region denotes the 
unitary point where l/(fcpa) = 0. 



that \/{kpa) increases with decreasing density and becomes 
positive after p/po ~ 0.06. This turning point is called the 
unitary limit, which we indicate by a red vertical line in the 
figure. We will discuss unitary limit in next subsection. 
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FIG. 3: (Color online) The critical temperature T c scaled by the 
Fermi energy sf as a function of density. Also shown is the BCS 
prediction (dashed line). 



The numerical result for the critical temperature is shown 
in Fig. 3. The dashed line shows the critical temperature 
over Fermi energy ratio given by BCS theory, which blows 
up quickly when density goes down in the crossover and 
BEC regions {T® cs is almost ten times larger than T c at 
density p = O.OOOlpo)- Physically, the ratio T C BCS /e F as 
well as At=o je f measures the strength of the attraction be- 
tween neutron and proton. However, due to the lack of pair- 
ing fluctuation effect, in crossover and BEC regions the BCS 
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theory does not give correct critical temperature for super- 
fiuid/normal transition which is mainly determined by the 
bosonic degree of freedom in these regions. The solid line is 
for T c obtained from Eq. (3.30). We can see that the evolution 
of T c is smooth and the superfiuid phase transition is second 
order in the whole density region. Also, it can be seen that T c 
is not a monotonous function of density: There is a local max- 
imum in T c curve which is roughly located around the unitary 
point. One should notice that a similar local maximum also 
appears in the famous Nozieres— Schmitt-Rink approach for 
T c [3]. At low density limit, all the nucleons participate into 
the deuterons which are long lived at temperature lower than 
T c , the system is essentially a deuteron gas and the superfiuid 
is totally due to the BEC of deuterons. In this case, we have, 



at T c , 2ZA pg (T = T t 
T c , we get, 



„uncondcnse 
Pb 



d = p/2. Solving out 



T c = 



2tt 



4C(3/2) 



2/3 



(4.6) 



This is just the BEC transition temperature for boson of mass 
mb- Adopting that mt ss 2m, we arrived at the well-known 
result, 



0.218£f, 



(4.7) 



which coincides well with our numerical result. 

It should be stressed that the BCS critical temperature 
yBCS was f oun( j t b e a g OOC [ approximation for the pair 
dissociation temperature T* (above which the pairs are essen- 
tially dissociated by thermal motion of the participators) [36]. 
So it is a pseudogap dominated region in between T c and 

rpBCS 

Corresponding to the evolution of the critical temperature, 
it is indicative to see how the pseudogap evolves. In Fig. 4, 
we plot the zero temperature excitation gap A(T = 0) as well 
as the pseudogap A pg (T = T c ) at T c , It can be seen that 
at low density p < O.Olpcii A pg (T = T c ) is roughly equal to 
A(T = 0) (but they do not completely coincide) reflecting the 
strong coupled nature in this case; however, at high density, 
A pg (T = T c ) is much smaller than A(T = 0), indicating that 
the pairing fluctuation is not essential there and BCS theory 
can work well. In the following subsections, we will focus 
on the crossover and BEC regions where BCS theory is not 
applicable at finite temperature. 



B. Unitary matter 

As shown in last subsection, the np scattering length kpa 
is infinite at the unitary point p « 0.06/9o- We call the nuclear 
matter at this point a unitary matter. It is interesting because 
it exhibits universal behaviors, i.e., the physical properties of 
unitary matter are independent of the details of the interac- 
tions [48]. Hence, the unitary nuclear matter behaves just 
like unitary cold atomic Fermi gas which has been realized 
in laboratory through the Feshbach resonance. For unitary 
matter, the unique characteristic scale is given by the Fermi 
momentum kp, so we have pr=o = C £ f, ^t=o = J£f, 




FIG. 4: (Color online) The excitation gap A at zero temperature and 
the pseudogap at T c as functions of density. 



T c = asF, A pg (T = T c ) = Xef, etc, with £, 7, a, A, etc, 
being universal constants. In our GqG scheme, the univer- 
sal coefficients are given by £ sa 0.59, 7 rj 0.64, a w 0.26 
and A « 0.53. For comparison, we would like to list the 
values obtained by Monte-Carlo techniques: £ w 0.42 [49], 
7 « 0.50 [49], a w 0.157 [50] or 0.25 [51]. Our results are 
larger than Monte-Carlo values. It is easy to show that the 
energy per particle in unitary matter is E/N = ((E /N)f lee 
where (E/N)f rcc is energy per particle for free fermion gas. 
Moreover, the equation of state of unitary matter is the same 
as free fermion gas, e = 3P/2, with P being the pressure. 
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FIG. 5: (Color online) The superfiuid gap A s f (dashed line), pseu- 
dogap A pg (dotted line) and the total gap A (solid line) as functions 
of the temperature at unitary point p — O.OGpo. The BCS result 
(dot-dashed line) above T c is also shown. 



In Fig. 5, we plot A s f (in units of Ep, the same below), A pg , 
and A as functions of temperature for unitary matter. As a 
comparison, we also plot the BCS result Abcs above T c . As 
we can see from the figure, with decreasing temperature be- 
low T c , Apg(T) is a monotonically decreasing function from 
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its maximum value at T c and it essentially vanishes at T = 
roughly according to A pg (T) cx T 3 / 4 [see Eq. (3.29)], while 
A s f (T) and A(T) both increase monotonically and become 
coincident at T = 0. Such kinds of temperature dependence 
reflect the fact that the pseudogap is due to the thermally ex- 
cited pairs: When T grows higher and higher, more and more 
pairs are excited from the condensate and at T c all the con- 
densed pairs are thermally excited; after that the thermal mo- 
tion of nucleons begins to dissociate the pairs and hence A pg 
(more exactly, -ZA pg ) begins decreasing above T c . In addi- 
tion, the critical temperature T c is smaller than the BCS pre- 
diction which shows the fact that the pairing fluctuation tends 
to destroy the order of the system. Although the physical pic- 
ture is clear, our formalism cannot be applied to very high 
temperature, where the effects of finite life-time of the pairs 
become significant which are not included in our formalism. 
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FIG. 6: (Color online) The entropy density and the specific heat of 
unitary matter as functions of temperature. Unlike the weak-coupling 
BCS case, there is no jump for the specific heat at T c . 



One should note that A and its derivative dA / dT are con- 
tinuous at T c . This is very different from the BCS case, 
where dA^cs/dT is discontinuous at T C BCS . Such differ- 
ence could be reflected in thermodynamic quantities such as 
the specific heat. In Fig. 6, we illustrate the entropy den- 
sity and the specific heat cy for unitary matter. We compute 
cy through cy = Tds/dT which involves the derivatives 
dp/dT, dA/dT, dm b /dT, and dp b /dT, so it is a nontriv- 
ial calculation. As is well-known, for the weak-coupling BCS 
case, the specific heat has a jump Acy cx dA 2 /dT at T C BCS 
which reflects the sudden opening of the excitation gap. For 
unitary matter, however, we found a continuous cy at T c . This 
continuity of cy reflects the previous existence of the excita- 
tion gap above T c due to the pairing fluctuation. It may serve 
as an experimentally accessible signal for the existence of the 
pseudogap in the normal phase. 



C. Deuteron gas 

When density is very low, say p < 0.002/?o from Fig. 2, 
the symmetric nuclear matter is effectively a Bose gas of 
deuteron. We in this subsection study the properties of 
deuteron gas based on our GoG formalism. First, we observed 
from Fig. 2 and Fig. 3 that when p — » 0, — \i T c . In this 
case, A, p and Z are almost temperature independent below 
T c reflecting the strong np attraction. Then for T < T c , the 
governing equations become (expanding in powers of a 3 p and 
A 2 /ii 2 ) [38], 




These relations give how A(p), A pg (p), p(p), and rrib(p) 
evolve with p at T < T c in the deep BEC region. 

Next, let us study the temperature dependence of these char- 
acteristic quantities. To specify the problem, we fix the den- 
sity as p = O.OOlpo- By solving the coupled equations (3.30), 
we get the transition temperature T c ~ 0.22 £p « 0.08 MeV. 
In Fig. 7, we show the superfluid gap A s f , the pseudogap 
A pg , the total excitation gap A as functions of temperature. 
Due to the stronger attraction, unlike for the unitary matter, 
now A(T) is almost a constant below T c , and at T = 0, A/ep 
is even larger than 1. But near zero temperature, A pg (T) still 
behaves as A pg oc T 3 / 4 as shown in Eq. (4.9). 

In Fig. 8, we give the temperature dependence of nucleon 
chemical potential p(T) and effective deuteron chemical po- 
tential Pb{T). Below T c , pt is zero meaning that BEC su- 
perfluid is formed. Above T c , both p and p b decrease, cor- 
responding to the thermal dissociating effect. Just above 
T c , simple calculation leads to that p(T) — p(T c ) oc pt cx 
(T — T c ) 2 . Since pt is related to p and the deuteron binding 
energy E\, through p b = 2/i + E b , we get that the binding en- 
ergy Eb at zero temperature is roughly 0.6ef at p = O.OOlpo 
from Fig. 8. The effective deuteron mass m^T) is shown 
in Fig. 9. It is seen that at low temperature mt is almost a 
constant, but it drops when temperature becomes higher. We 
note here that mj can be regarded as the medium renormalized 
deuteron mass only in deep BEC region [52, 53]. It is renor- 
malized because it contains indirectly the deuteron-deuteron 
interaction through the nucleon-deuteron coupling in the nu- 
cleon self-energy. Hence, this effective mass is not equal to 
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FIG. 7: (Color online) The superfluid gap A s f, pseudogap A pg 
and the total gap A as functions of the temperature at density p = 
O.OOlpo. The BCS result is also shown above Tc- 



FIG. 9: (Color online) The effective deuteron mass parameter as a 
function of T at density p = 0.001 po. 



pairs. The former contribution is just the BCS theory result, 



2m — Ef,, as one may intuitively expect. It is actually a param- 
eter measuring the effective size of the non-condensed pairs, 
hence it is also defined in intermediate coupling and even 
weak coupling regions. The drop of m& at high temperature 
simply indicates that the effective size of the non-condensed 
pair is enlarged by the thermal motion of participate nucleons. 
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FIG. 8: (Color online) The nucleon chemical potential and the ef- 
fective deuteron chemical potential as functions of T at density 
p = 0.001/00. 



('^) 1/2 e- Ao/T , T«T C 



(4.10) 



with Ao = A<r=o- The latter one is dominated by 
TdAp g /dT for other quantities are almost independent of T 
[see Eq. (4.9)], hence 



T 3/2 j, < j, 



(4.11) 



which dominates cy at low temperature. At the phase transi- 
tion point, similarly with unitary matter, due to the continuity 
of the temperature derivative of excitation gap, cy does not get 
a discontinuity, but a A-type behavior. Now, both the low tem- 
perature and the critical behaviors of cy are quite similar with 
the situation found in ideal BEC superfluid. It indicates that 
the symmetric nuclear matter at very low density is a nearly 
ideal deuteron gas. 

Figures. 6 and 10 inspire us that the specific heat jump at 
T c may serve as a possible thermodynamic signal for BCS- 
BEC crossover. We hence draw in Fig. 1 1 the specific heat 
jump Acy = cy (T c - 0+) - cy (T c + 0+) at T c over density 
ratio as a function of density. As density decreases, Acy / p 
monotonously decreases and vanishes when density is smaller 
than 0.06po which is just the unitary point. The physical rea- 
son for such kind of behavior is clear: as density decreases the 
system becomes more and more bosonic and the finite jump 
of the specific heat at T c which is a typical BCS feature gets 
suppressed. 



Finally, we depict the entropy density and the specific heat 
for deuteron gas in Fig. 10. One should note that at low tem- 
perature, the behavior of cy is very different from the predic- 
tion of BCS theory: It shows T 3//2 dependence at low T rather 
than an exponential suppression. Actually, since the conden- 
sate does not contribute to entropy, cy at T < T c contains 
contributions from quasi-nucleons and from thermal excited 



V. SUMMARY AND DISCUSSION 

As is well known, BCS theory is only applicable to weak 
coupling system or at zero temperature since it does not 
contain any pairing fluctuation effects. To study the BCS- 
BEC crossover at finite temperature, it is necessary to go be- 
yond BCS description. We, in this article, studied the ef- 
fects of pseudogap due to pairing fluctuation on the BCS-BEC 
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FIG. 10: (Color online) The entropy density and the specific heat 
of symmetric nuclear matter as functions of temperature at density 
p = O.OOlpo. 
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FIG. 11: (Color online) The specific heat jump at T c over density 
ratio Acv / p as a function of density. 



crossover problem in symmetric nuclear matter. For this pur- 
pose, we adopted a T-matrix method based on a GqG approx- 
imation for the pair susceptibility. This method is a natural ex- 
tension of BCS theory and has been widely used in theoretical 
studies of high T c superconductor and BCS-BEC crossover 
problems in cold fermion atoms. 

The pseudogap is determined by the density of thermally 
excited rip pairs and vanishes at zero temperature. We found 



that its effects are substantial for intermediate and strong cou- 
pling regions (corresponding to intermediate and low-density 
regions) in the BCS-BEC crossover when temperature is not 
zero. At high-density region, the pseudogap is essentially 
small, and the GqG theory recovers the BCS theory. Taking 
into account the pseudogap effects, we calculated the critical 
temperature for superfluid phase transition shown in Fig. 3. 
At high density, T c follows the BCS result, but at interme- 
diate and low densities, it deviates from the BCS prediction 
remarkably. At dilute limit, T c coincides with the BEC transi- 
tion temperature of dilute deuterons. 

The pseudogap persists in both T > T c and T < T c re- 
gions. We investigated how the pseudogap affects the prop- 
erties of unitary matter and dilute deuteron matter. At inter- 
mediate and low densities, the significant result is that due to 
the pseudogap, the specific heat is continuous at T c . At low 
density, cy oc T 3 / 2 at low temperature and has a continuous 
A-type behavior at T c just like an ideal BEC superfluid. The 
qualitative change of the temperature dependence of specific 
heat from high density to low density also indicates the BCS- 
BEC crossover. Moreover, as density decreases, the jump of 
specific heat at T c decreases and eventually disappears when 
p < 0.06po as shown in Fig. 11. This may serve as a thermo- 
dynamic signal for BCS-BEC crossover. 

We stress that the GqG approximation is, in principle, not 
applicable at temperature much higher than T c , since then the 
thermal dissociation effect will result in a finite width of the 
pairs which we did not take into account. Besides, symmetric 
nuclear matter may favor the formation of a cluster [54], but 
we did not consider this situation in present article. Finally, 
since most of the nuclear systems in nature do not contain 
equal numbers of neutrons and protons, it will be significant to 
extend present formalism to asymmetric nuclear matter [14]. 
We leave all these challenges for future studies. 
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